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KREEP materials were thought to be last crystallized at the lunar crust and mantle boundary. Impact 
cratering and volcanism are mainly responsible for their distributions on the lunar surface. Therefore, 
observation of global KREEP materials and investigation of distributions in the areas of large basins are of 
critical importance to understand the geologic history of the Moon. Here we report the new global 
potassium distribution on the Moon detected by Chang'E-2 Gamma-ray Spectrometer. We found that our 
new measurements are in general agreement with previous observation. A new finding and an important 
difference is that relatively higher K abundances in the Mare Crisium and Mare Orientale than their 
surrounding rims were detected for the first time. In light of our observations in these two areas, we propose 
that Crisium and Orientale basin-forming impact events may have penetrated to the lower crust and 
excavate the deeper materials to the lunar surface. 



I n the Lunar Magma Ocean (LMO) hypotheses 1 , KREEP materials are proposed to lie between the mafic mantle 
I and plagioclase crust in the process of compositional crystallization. Their presence on the lunar surface is the 
I end product of a series of evolutionary processes that have changed the surface features of the Moon 2 . In light 
of LMO hypotheses and the general impact cratering theory 3 , impact events penetrating through the crust should 
excavate the KREEP materials on the lunar surface. Therefore, observation of the global KREEP material on the 
Moon and the investigation of KREEP material in the areas of large basins are of critical importance to understand 
the formation of such a magma ocean and constrain the depth of the huge impact events. 

In the past few years, several spacecraft with improved compositional remote sensing instruments were sent 
into the orbits around the Moon and large quantities of valuable data have been obtained. Chang'E-2 spacecraft, 
the second lunar probe of China, was launched on 1 October 2010, and inserted into a circular orbit of 100 km 
altitude on 15 October 2010. It carried a suite of remote sensing and plasma instruments. Here we report the 
results from the gamma ray spectrometer during its 178-days measurements and their new implications on the 
impact depth of the Crisium and of Orientale basin, respectively. 

CE-2 GRS is the first gamma ray spectrometer that used large LaBr 3 crystal as the main detector to measure 
gamma rays emitted from the Moon within the energy region from 0.3-MeV to 9-MeV 4 ' 5 . The energy resolution of 
the flight model is — 5% full width at half maximum at 662 keV, which is about two times better than LPGRS 6 . 
Since CE-2 GRS was designed to be mounted inside the spacecraft, a big well-shape CsI(Tl) crystal was employed 
in anticoincidence with the main crystal to suppress the Compton effects and the background gamma rays from 
the spacecraft materials generated by Galactic Cosmic Rays (GCRs). Gamma rays emitted from the lunar surface 
are mainly from natural radioactive elemental (e.g., U, Th, K) decay effects and major elemental (e.g., Si, O, Fe, Ti, 
Mg, Al, Ca) excitation from interactions with incident GCRs 7 . Gamma rays entering the main crystal were 
recorded by the CE-2 GRS with a cadence of 3 second in the 512-channel spectrum. 

More than 25 full energy peaks could be identified in CE-2 GRS spectra corresponding to different elements on 
the lunar surface (Fig. 1). But owing to the limitation of intrinsic energy resolution of the scintillation detector, 
many spectral peaks contain contributions from two or more elements. The strongest peak around 1.46 MeV is a 
combination of lines resulting mainly from the K line at 1.46 MeV and LaBr 3 self- activity because of the 138 La 
decay (1.436 MeV + 32 keV) 8 that contributes about 85% of the corresponding gamma-ray intensity. 
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Figure 1 | Chang'E-2 gamma-ray spectrum plotted as counts per second (cps) of the whole Moon versus energy for its 178-days measurements with 
100 km altitude. The major peaks are labeled. The spectrum below 2.7 MeV are mainly contributed by the self-activities from LaBr 3 crystal: natural 138 La 
is primarily responsible for the spectrum below 1.7 MeV while 227 Ac decay chain generates alpha emitters between 1.7 and 2.7 MeV, as label 3. K line 
( 1 .46 MeV) is partly contributed to the sum peak as Label 2 with the 138 La decay to 138 Ba ( 1 .436 MeV + 32 keV) of the LaBr 3 crystal; label 1 is the sum peak 
of elements from the Moon, and crystal beta decay of 138 La (789 keV) to 138 Ce. The inset figures are the fitting results of the energy region around 
1.46 MeV. (A) The orbital spectrum and its fitting result. (B) The cruise spectrum and its fitting result. 



Corrections of gain, solid angle, and cosmic ray variation have 
been carried out to CE-2 GRS measurements. The spatial resolution 
of CE-2 GRS is about 150 km in its 100 km working orbit, but in this 
work, all the corrected spectra are binned into equal-area pixels about 
60 km X 60 km on the Moon. In our calculation we assumed gamma 
rays emitted from LaBr 3 self-activity decay are constant in the whole 
mission, the spectra measured in the cruise period could therefore be 
used as the background. The net 40 K gamma- rays emanating from 
the Moon in each pixel are obtained by calculating the peak fitting 
difference between binned spectra and the cruise spectrum (see inset 
of Fig. 1). Then all the net 40 K count rates from 11306 pixels were 
compared to count rates in the detector derived from the calculated 
surfaced gamma-ray fluxes to determine the abundance 7,910 . 

Results 

The K distribution from CE-2 GRS (Fig. 2) is in general agreement 
with Lunar Prospector 91112 , Kaguya 10 , and Chang'E-1 measure- 
ments 13 . The global average abundance is 620 ±615 ppm with the 
upper limit of 3240 ppm, which is lower than the LP value of 755 ± 
740 ppm with the upper limit of 4360 ppm 9 . The Procellarum 
KREEP Terrane (PKT) 14 has an average value of 1768 ± 530 ppm, 
comparing to 2004 ± 708 ppm determined by LPGRS. In the South 
Pole-Aitken (SPA) basin, the abundance ranges from 500 ppm to 
1350 ppm with an average value of 615 ± 235 ppm, which is also 
lower than the LP measurements of 765 ± 320 ppm. Different from 
the PKT and SPA basin, lunar highlands is depleted in the K with the 
abundance being less than 500 ppm according to the CE-2 GRS 
observations. This is consistent with LP measurements 9 . 

Mare basins on the nearside, such as Serenitatis, Tranquillitatis, 
Fecunditatis, and Crisium, have lower K abundance than mare 
regions. But after careful examination of the CE-2 GRS data, we 
found that the abundance within the Crisium impact basin is clearly 
higher than its surrounding highlands (Fig. 3a-b). This is inconsist- 
ent with previous LP results (Fig. 3c), which indicates that there is no 



distinct difference between the K abundance of Mare Crisium and 
that of its surrounding 11 . 

Discussion 

The Crisium basin is a Nectarian-age multi-ring impact structure on 
the eastern nearside of the Moon. Its interior is covered by mare 
basalt from a protracted period of volcanism over 0.8 Ga 15 . The 
abundance in its inner mare basalt ranges from 430 ppm to 
840 ppm with an average value of 600 ±100 ppm. On the other 
hand, the K abundance in the surrounding highlands is lower than 
480 ppm with an average value of 290 ± 105 ppm. Because the 
southwest highlands is thought to be covered with ancient mare 
deposits 16 , not with ejecta from the Crisium basin, the abundance 
in its northeast highlands thus represents the true value of the ejecta 
for estimating the average value here (Fig. 3b). The highest K abund- 
ance in the Picard crater (14.6°N, 54.7°E) indicates the subjacent 
material has K abundance about 3 times (900 ppm) higher than 
highlands deposits surrounding the Crisium basin. 

There is a long standing question concerning the origin of material 
excavated by the Crisium impact event 16,17 . Some studies indicated 
that the Crisium impact event only excavated the upper or middle 
crustal material with a depth of about 40 km 16 while others preferred 
that the impact event might have excavated mantle material 17 . 
According to recent observations of olivine sites around the 
Crisium basin 17 , the scenario that impact event only excavate upper 
or middle crustal material is probably excluded because olivine is a 
significant component associated with lunar mantle or Mg-rich plu- 
tons in lower crust 18 . Furthermore, if the Crisium impact event did 
penetrate through the crust and KREEP layer, exposing the upper 
mantle, KREEP (e.g., K here) and mafic materials should be found at 
its rim in light of the general impact cratering theory 3 , deep material 
will be spread around the near rim. However, higher K signal is not 
observed in Crisium rim and its basin terrain by CE-2 GRS, FeO is 
also found as lower as 2-4 wt.% in these areas 19 . 
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Figure 2 | The global K abundance (unit: ppm) of the Moon from CE-2 GRS 178-days measurements. All the data are binned into quasi-equals 2° X 2° 
pixels in the equator and displayed in a cylindrical projection. 



On the basis of our observation and previous results 161719 , we 
prefer the scenario that the Crisium impact event only penetrated 
to the plutonic lithosphere resulting from the magmatic intrusion 
into the lower crust. The magmatic intrusion, which is enriched in 
the PKT, is the hybrid source with low- density Mg-rich olivine at the 
base crust resulting from the overturn of the early mantle cumulate, 
KREEP, and crustal anorthosite 20 . While the Crisium impact event 
penetrated through the crust, to the upper boundary of the magmatic 
intrusion, the deepest materials would be ejected to the surrounding 
rim area. The redistribution of excavated material resulted in the 
lower K and FeO in the highlands 19 bounding the basin and is in 
agreement with the observations of olivine exposure on the basin 
rim 17 . The constraint of impact to the upper boundary of the mag- 
matic intrusion is consistent with the higher K abundance in sub- 
jacent highlands material covered by mare basalt as detected in the 
Picard crater, and is also consistent with the non-zero crustal thick- 
ness in the Crisium basin 21 . The big basin formation impact exca- 
vated tremendous volume material and formed a thinner crust with 
fracturing underlying lithosphere 22 in the Crisium area. Owing to the 
density difference between lower crust and partly melted magma, the 
buoyancy contrast 23 enabled the eruption of basalts at the Crisium 
basin. The partly melted magma is from the Moon's deep interior 
and produced by slowly heating up from the decay of radioactive 
elements. This supports the observations of higher K abundance 
(about 600 ppm) and very higher FeO contents 19 (up to 20 wt.% 
for areas of the central areas) in the mare basalts covered in the 
Crisium basin. Located near the PKT region and consistent K abund- 
ance with mare basalts in the PKT boundary, high-heat production 
by radioactive decay is suggested to be responsible to the long dura- 
tion of volcanic extrusions at the Crisium basin 16 . 

Similar characteristics are also found in the Orientale basin 
(Fig. 4a-d), the youngest multi-ring basin on the Moon 24 . CE-2 
GRS shows that the central areas (5-6 pixels in diameter with each 
pixel of ~ 60 km) of the Orientale basin corresponding to the Mare 
Orientale, the largest mare deposit (~ 320 km in diameter) in the 
basin, has obviously high K abundance with the average of 360 ± 



115 ppm (Fig. 4b). This value is about a factor of two higher than the 
surrounding areas, including the Montes Rook Formation (MRF) 
occurring mostly between the edge of the Inner Depression (ID) 
and Cordillera Ring (CR), and Hevelius Formation (HF) lying out- 
side the CR. The MRF, consisting of smooth plains and corrugated 
deposits, has K abundance of 115 ± 95 ppm (Fig. 4b), which is only 
one-third of that in the Mare Orientale. Most areas of HF, which 
represent the radial ejecta deposits from the Orientale basin, have an 
average abundance a little higher (about 80 ppm) than the deposits 
in MRF, except two areas located at the southwest (K-l) and north- 
west (K-2) portion of the basin (Fig. 4d). These two areas, related to 
the well-preserved ejecta units in MRF (part of K-l) and HF 25 , have 
abundance of 270 ± 80 ppm (K-l) and 255 ± 90 ppm (K-2), 
respectively. Because turbulent impact modification after Orientale 
impact event has less influence on these two areas 26 , the K abun- 
dances in these two areas therefore represent the original value of 
excavated ejecta from the basin-forming impact. 

Assuming the Out Rook Ring is the closest approximation to 
the basin excavation cavity, CR is believed to have been formed by 
the inward collapse at the edge of structural uplift along the basin 
of the displaced zone 25 . MRF therefore represents the modification 
of rim ejecta and the mixture with the upper crust, which could be 
consistent with the similar K abundance in HF and MRF. The 
relative higher K abundance (about 80 ppm) in HF implies that 
the deeper material excavated by the impact event has relatively 
higher K abundance. The LMO scenario proposes that the KREEP 
layer lies between the anorthositic lunar crust and mafic cumu- 
lated mantle in the long-time compositional differentiation 1 . 
Therefore, in contrast to previous works 27 which suggested that 
Orientale impact event excavated only to the upper half of the 
lunar crust in its target region with the depth of 40-50 km, we 
prefer that the Orientale impact event penetrated into a source of 
radioactive- rich material at the lower crust. This scenario supports 
the lower FeO (2-4 wt.%) 19 in the surrounding of the Orientale 
basin and is also consistent with the estimated crustal thickness in 
this area 21 . This scenario is also consistent with the enrichment of 
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Figure 3 | Surface topography and K distributions of the 550-km-diameter Crisium Basin. All the maps are displayed in a Lambert azimuthal equal-area 
projection superimposed on the shaded lunar surface relief. Two circles represent the proposed rims with diameter of 300 km and 700 km. NE - northeast 
highlands of Crisium basin; SW - southwest of Crisium basin that covered by ancient mare basalt. Picard - Picard crater with 30 km in diameter. (A) the 
topographic elevation of the Crisium basin from Lunar Orbiter Laser Altimeter 32 . (B) K abundance of the Crisium Basin measured by CE-2 GRS with data 
binned into quasi-equal 2° X 2° pixels. (C) K abundance of the Crisium basin measured by LPGRS 9 . 



low-Ca pyroxene in the central peaks of the Maunder crater, a 
post-Orientale crater with 55 km diameter in the interior depres- 
sion of the Orientale basin, which is thought to be sampling the 
noritic lower crust 26 . 

The internal accumulation of radiogenic heat sources causing par- 
tial melting of the mantle is the basic postulation for the formation of 
mare basalt 28 . The details of the generation, ascent and eruption of 
the mare basalts in the Orientale basin need further investigations, 
but the similar K abundance in the Mare Orientale to the nearside 



mare basalts, especially to the boundary of PKT, suggests that mare 
deposits within Orientale basin are similar in composition to near- 
side mare basalts 27 . Furthermore, two areas with higher K abundance 
in the Mare Orientale definitely correspond to the eastern and north- 
western region of the inner depression covered by mare basalts, 
which have almost the identical age, ~ 3.2-3.3 Ga 29 ' 30 . The same 
concentration of FeO, Ti0 2 19 , and K in these two areas supports 
the interpretation that the covered basalts in the center of Orien- 
tale basin are from single volcanic events 29,31 . 
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Figure 4 | Orientale Basin and the K distribution from CE-2 GRS and LPGRS in cylindrical equidistant projection superimposed on the shaded lunar 
surface relief. (A) the topographic elevation of the Orientale basin from Lunar Orbiter Laser Altimeter 32 with the black circles represent the rings of the 
Orientale basin (CR, ORR), but the Inner Rook Ring is not shown here; CR - Cordillera Ring; ORR - Out Rook Ring; ID - Inner Depression; MRF - 
Montes Rook Formation; HF - Hevelius Formation; M-R - Mendel-Rydberg basin. (B) K abundance in the Orientale Basin measured by CE-2 GRS with 
data binned into quasi-equal 2° X 2° pixels. (C) K distribution in the Orientale basin measured by LPGRS 9 . (D) the smoothed K distribution from CE-2 
GRS measurements; K-l/2 represent two areas with relatively higher K Abundances outside the CR. The yellow lines represent the contours of K 
abundance in K-l/2 area. 



Methods 

All raw data received from the spacecraft were converted to time series gamma ray 
spectra that were used to derive elemental maps. However, in order to ensure the 
correctness of spectral analysis, mapping, and abundance calculations, data selection 
was carried out at first. We eliminated the data which have bad sync words, bad check 
sums, and off-limit values. The measurements during solar energetic particle event 
were also excluded. In this work, the energetic proton data taken from the Earth 
orbiting GOES satellites 33 were used to identify the solar energetic particle events. 

Three major corrections to the time series data were made for the variations of 
gain, galactic cosmic ray, and solid angle before analyzing the spectra. To monitor 



the gain changes, we calculated the peak positions of two strong peaks 
(0.511 MeV, 7.60 MeV) for every 15-min and every 3-hour interval. The time 
series data were then corrected according to the gain and offset values for each 
3 -hour and 15-min interval 34 . The time- varying galactic cosmic ray flux was 
monitored by measuring the counting rate of oxygen gamma ray line (6.13 MeV) 
for all latitudes poleward of + 85°. In order to apply the galactic cosmic ray 
correction, the 6.13 MeV counting rate were smoothed and normalized to the 
counting rate seen at the beginning of the orbital mission. All spectra were then 
divided by this normalized flux for the correction of the galactic cosmic ray 
variations 34 . In addition, all spectra were normalized to a spacecraft height of 
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100 km to normalize the variations of solid angle seen by the detector due to 
height variation from 85 to 125 km. 

All corrected spectra were binned into approximately equal-area pixels about 
60 km X 60 km by using their latitude and longitude positions at the times they were 
taken. Gamma-ray spectra were accumulated in each pixel to present the spectrum 
measured from such area. To calculate the K counting rate in each pixel, the inte- 
grated area beneath K peak should be determined. In this work, peak fitting method 
(Gaussian shaped, linear background line) was used to determine K peak area in each 
pixel. Since the peak around 1 .46 MeV is a combination of lines resulting mainly from 
the K line at 1.46 MeV and LaBr 3 self- activity because of the 138 La decay (1.436 MeV 
+ 32 keV) 8 , the contribution of LaBr 3 self-activity should be eliminated from each 
calculation. We assumed that gamma rays emitted from LaBr 3 self- activity decay were 
constant in the whole mission. The peak area estimated from the accumulated 
spectrum in the cruise period was used to eliminate the contribution of LaBr 3 self- 
activity. The net 40 K gamma-rays emanating from the Moon in each pixel were 
obtained by calculating the area difference between binned spectra and the cruise 
spectrum. 
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